Dynamics of Feshbach Molecules in an Ultracold Three-Component Mixture 
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We present investigations of the formation rate and collisional stability of lithium Feshbach 
molecules in an ultracold three-component mixture composed of two resonantly interacting fermionic 
"^Li spin states and bosonic ^^*Yb. We observe long molecule lifetimes (> 100 ms) even in the pres- 
ence of a large ytterbium bath and extract reaction rate coefficients of the system. We find good 
collisional stability of the mixture in the unitary regime, opening new possibilities for studies and 
probes of strongly interacting quantum gases in contact with a bath species. 



Magnetic Feshbach resonances allow precise control of 
collisional properties, making them a key tool in ultra- 
cold atom systems. They have been extensively used to 
study ultracold molecules, as well as few- and many-body 
physics Two-component Fermi gases near a Fesh- 
bach resonance provide excellent opportunities to study 
strongly interacting quantum systems 0- This is possi- 
ble due to the remarkable collisional stability of the atom- 
molecule mixture on the positive scattering length side 
of the resonance [1, 01 , attributed largely to Fermi statis- 
tics Q. Extending the system to three-component 
mixtures in which only two are resonantly interacting 
offers the exciting possibility of modifying or probing 
pairing dynamics by selective control of the third compo- 
nent. A third-component may also be used as a coolant 
bath for exothermic molecule-formation processes, pro- 
vided that inelastic processes with the bath are negligi- 
ble. In the context of many-body physics, a third non- 
resonant component can be useful as a microscopic probe 
of superfluid properties 0, |1] ; as a stable bath for stud- 
ies of non-equilibrium phenomena 0], or as an accurate 
thermometer of deeply degenerate fermions [T^] ■ 

Collisional stability of Feshbach molecules in the ab- 
sence of Fermi statistics becomes a crucial question for 
multi-component mixtures 0, [lil, [HI • A recent theoret- 
ical analysis of such mixtures suggests a possibility for 
enhanced molecule formation rates with good collisional 
stability [ll|. Enhanced atom loss has been observed 
near a ^Li p-wave resonance in the presence of a ^^Rb 
bath [31, while a small sample of the probe species 
has been found to be stable within a larger strongly in- 
teracting ^Li sample [3- 

In this paper we investigate a mixture composed of 
two resonantly interacting spin states of fermionic ^Li im- 
mersed in a large sample of bosonic ^^^Yb atoms. While 
the Li interstate interactions are arbitrarily tunabl e by 
means of an s-wave Feshbach resonance at 834 G [14| . 
the interspecies interactions between Li and Yb are con- 
stant and small We present the first observations 
of formation and evolution of Feshbach molecules in a 
bath of a second atomic species. In the unitary regime, 
we observe good collisional stability of the mixture with 
elastic interactions dominating over inelastic losses. We 
extract the reaction rate constants from a classical rate 
equations model of the system. 

Our experimental procedure has been described in ear- 



lier work [11]. Briefly 3 x 10** atoms of ^'^"Yb in the 
^5*0 state and up to 4 x 10^ atoms of ^Li , distributed 
equally between the two ^6*1/2, F = \ states (denoted 
Li|l), Li|2)), are loaded from magneto-optical traps into 
a crossed-beam optical dipole trap. We then perform 
forced evaporative cooling on Yb to the final trap depth 

t'Yb(Li) = 15(55) ^KjWith mean trap frequency a}Yb(Li) — 
2tt X 0.30(2.4) kHz [17[, during which Li is cooled sympa- 
thetically by Yb. Following evaporation, the mixture is 
held at constant trap depth to allow inter-species ther- 
malization. With a time constant of 1 s the system ac- 
quires a common temperature Tyb = Tli = 2 fiK with 
atom number iVYb(Li) = 2 x 10^(3 x lO"). This corre- 
sponds to Ti^i/Tp ~ 0.4, and TVb/Tc ~ 2.5, where Tp is 
the Li Fermi temperature and Tc is the Yb Bose-Einstein 
condensation temperature [ist . 

After this initial preparation, we ramp up the mag- 
netic field to a desired value and observe the system af- 
ter a variable hold time. For fields in the vicinity of the 
Feshbach resonance, there is a field-dependent number 
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Figure 1: Li atom loss spectroscopy in the presence (filled 
circles) and absence (open squares) of an Yb bath near the 
^Li 834 G Feshbach resonance (inset). We plot the number 
of Li atoms after 500 ms of evolution normalized to that at 
10 ms. The thick dashed line indicates the resonance center 
and the thin dashed line indicates the magnetic field at which 
~ ksTi^i for the initial conditions. 
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loss and heating for the Li cloud during the 20 ms ramp 
time, resulting in Tl; rising to as high as 4.5 /iK. At this 
point, the density- weighted average density (riYb(Li)) is 
2.6(0.35) X lO^^cm"^. For interrogation in the absence 
of the bath, Yb is removed from the trap with a 1 ms 
light pulse resonant with the ^5*0 — !> ^Pi transition [l9l |. 
Atom number and temperature are monitored using ab- 
sorption imaging for both species after switching off the 
magnetic field. 

We first present our results on atom loss spectroscopy 
near the Fcshbach resonance (see Fig.[T]). The atom- loss 
maximum obtained in the absence of Yb has been ob- 
served previously [l^ and can be explained as a result 
of the formation and subsequent decay of shallow lithium 
Fcshbach dimers [sl. 13. [20l - l22j which form only on the pos- 
itive a side of the resonance. Here a denotes the Li|l)- 
Li|2} scattering length. In the presence of the Yb bath, 
the loss feature is shifted and broadened. We interpret 
the behavior of the mixture in terms of five chemical pro- 
cesses: 
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Forward process (I) corresponds to a three-body collision 
event which produces a shallow Fcshbach dimer (denoted 
Lij) accompanied by the release of the dimer binding 



energy 



Li denotes a Li atom in either of 
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the two spin states. Process (II) corresponds to two-body 
loss to a deeply bound dimer (denoted Li2) with binding 
energy en- Processes (III) and (IV) are similar to (I) 
and (II) with the spectator atom being Yb rather than 
Li [2j|. Process (V) corresponds to direct three-body 
loss to a deeply-bound molecule. Processes (II, IV, V) 
always result in particle loss from the trap since to ^ 
L^Li- Vibrational relaxation due to collisions between Li2 
Fcshbach molecules may contribute at the lowest fields, 
but has a ne gligibl e rate for the fields at which we perform 
our analysis [3, [22| • We have experimentally checked that 
direct three-body loss processes to deeply-bound states 
involving three Li atoms as well as those involving one Li 
atom and two Yb atoms are negligible for this work [24 1 . 
Three-body losses involving Yb atoms alone have a small 
effect (25} and are taken into account in our analysis. 

In the absence of Yb, only processes (I) and (II) con- 
tribute. If we neglect loss process (II), the atom- molecule 
mixture approaches an equilibrium, characterized by an 
equality of the forward and reverse rates and an equilib- 



rium molecule fraction jq^._^2N 
where Nm is the molecule number and (pi^i is the phase 
space density for each spin component in the ground 
state of the trap [20l - l22j . The timescale for achieving 
equilibrium depends on the three-body rate constant L3 
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Figure 2: Evolution of Li Feslibach molecule number at 709 
G without (a) and with (b) an Yb bath. The numbers are ob- 
tained by comparing Li atom numbers (insets) ramped across 
resonance (diamonds) or not (open squares) as described in 
the text. Lower inset also shows Yb number (filled squares). 
The curves are fits with a rate equations-based model. 



for process (I), which scales with the scattering length 
as a^, whereas rate constant L2 for process (II) scales 
as a~^-^ [2^, [23] . The shape of the loss spectrum can 
thus be qualitatively explained by noting that the dimer 
formation rate increases with magnetic field while equi- 
librium dimer fraction and molecule decay rate decrease. 
The large rate for process (I) at high fields close to reso- 
nance ensures equilibrium molecule fraction at all times. 
Broadly speaking, the rate-limiting step determining the 
system evolution is the molecule formation rate at low 
fields and decay rate at high fields. The trap depth 
also affects the loss spectrum shape, since it determines 
the magnetic field range over which the formed shallow 
dimers remain trapped. 

In the presence of Yb. the additional dimer formation 
(III), dimer decay (IV), and 3-body loss (V) processes 
contribute. The observed loss spectrum is broadened on 
the higher field side, suggesting that for our parameters, 
processes (IV) and/or (V) play an important role while 
process (III) does not. The rate constants L3, L'2, and 
L3, for processes (III), (IV) and (VL have theoretical 
scalings a^, a~^, and a^, respectively [lll,!!^. Overall, we 
see two regimes of behavior - a lossy one where molecule 
formation is energetically favored (es > /cbTli) and a 
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stable one closer to resonance (e^ < ^sTLi). The crite- 
rion eg = kgT separating these two regimes is equivalent 



to fca = 1 where 



2inLi 



ksTu, i-C-, the unitary criterion. 



In order to expand upon this qualitative picture, we 
study the time evolution of the three-component mix- 
ture at representative magnetic fields in the above two 
regimes. We are then able to extract quantitative in- 
formation for the above processes from a rate-equations 
model of the system. 

Fig. [2] shows the Li atom and molecule number evo- 
lution at 709 G (es — ks x 8.3 /iK) a field value where 
modifications due to the Yb bath are apparent in Fig.[T] 
The number of Feshbach molecules at a particular field is 
determined by using a procedure similar to earlier works 
|3l ■ After variable evolution time, we ramp the mag- 
netic field with a speed of 40G/ms either up to 950 G, 
which dissociates the molecules back into atoms that re- 
main in the trap, or to 506 G, which does not. We then 
rapidly switch off the magnetic field and image the atomic 
cloud. The molecule number is obtained from the num- 
ber difference in the two images (see insets in Fig. [2]). 

We sec that the presence of Yb alters the molecule 
decay rate while the formation rate is unchanged. The 
Feshbach molecules appear to coexist for a long time 
(> 100 ms) with the Yb bath, even in the absence of 
Pauli blocking We adapt the recent rate-equations 
analysis of Feshbach losses in a Fermi- Fermi mixture [22| 
to incorporate a third component, temperature evolu- 
tion, and trap inhomogeneity. Tli/Tp > 0.5 is satisfied 
throughout the measurement range, allowing a classical 
treatment of the Li cloud. We model the density evolu- 
tions due to processes (TV) using: 

rim = Rm+ R',n ^ L2nmnu - L'^nraUYh (1) 

riLi = -2i?,„ - 2i?J„ - L2n,nnu ~ 2Lfn'linYh (2) 
"Yb = -L2n„inYb ~ LfnliUYb. (3) 

Here , riLi and riyb sltc the densities of shallow dimers 
Li|, Li atoms and Yb atoms, respectively. Rm{R'm) — 
\Lz{L'^)n'l(ni^i{Yh)-qL'i{L'r^)nrani^\{Yh) is the net-rate for 
molecule production via process (I) ((III)). We determine 
q through the constraints on the molecule fraction at 
equilibrium {Rm{R'ra) = 0)- We obtain an upper bound 
for L3 by observations at large negative a (described be- 
low) which indicates a negligible effect for the data in 
Fig. [21 allowing us to set = for the analysis at 709 G. 

The time evolution of Tli and Tyb are modeled consid- 
ering the energy deposition from processes (I) and (III) as 
well as heating from the density-dependent loss processes 
(II) . (IV) and (V) • In addition, our model also takes 
into account the effects of evaporative cooling [2^ , inter- 
species thermalization [l5{ . one-body losses from back- 
ground gas coUisions, and Yb three-body losses [25l. [28j. 
The Li scattering length at 709 G is a = 1860 ao, ensur- 
ing rapid thermalization (< 1ms) in the lithium atom- 
Fcshbach molecule mixture [26} . This allows the assump- 
tion of equal temperature Tli for lithium atoms and Fesh- 
bach molecules. The heating from molecule-formation at 
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Figure 3: The evolution of temperature and number at 810 G 
for Li atomic cloud with Yb (filled circles) and without (empty 
circles) and also for Yb in the presence of Li (filled squares). 
The curves are fits with a rate equations-based model. 



709 G dominates over inter-species thermalization, main- 
taining ~ 4.5 //K and Tyb — 2^K, as observed in 
both experiment and model. 

The best-fit rate coefficients extracted from the atom 
data (shown in the insets) are L3 = (1.4 ± 0.3) x 
lO-^'^cmVs, £2 = (1.3 ±0.3) X IQ-^^cxa^/s, and L'^ = 
(2.3 ± 0.2) X IQ-^^ cm^/s. is consistent with 0. Ah 
reported uncertainties are statistical. The L3 value is 
consistent with that obtained in [^ , after accounting for 
the slight differences in experimental parameters. Using 
iji^-Yb) as a measure of the dimer decay rate, we get 
170 ms as the lifetime of a Li Feshbach molecule in the 
Yb bath. 

We now turn to the unitary regime, where we choose 
810 G (fca = +6, = fcs X 0.11 /iK) as our representa- 
tive field to study the mixture properties. It is difficult 
to reliably observe the molecule number using our ear- 
lier method in this regime, so we only monitor the atoms 
(see Fig. [3]). Starting with an inter-species temperature 
differential as before, we observe a fast drop in Tli in the 
presence of Yb and clear evidence of inter-species ther- 
malization. The Li number in the three-component mix- 
ture exhibits a long 1/e lifetime of 2 s, far larger than at 
709 G. However this is still an order of magnitude shorter 
than that obtained in the absence of Yb. The interpreta- 
tion of the decay is not straightforward as both two-body 
(process (IV)) and three-body (process (V)) inelastic loss 
can contribute [3, H^l ■ The large rate for process (I) in 
this regime ensures equilibrium molecule fraction at all 
times. By fitting to data taken at 935 G where ka — —2 
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and process (V) is expected to dominate inelastic loss, 
we obtain = (4.3 ± 0.3) x IQ-'^^cm^/s. This sets a 
lower bound for L3 at 810 G. We fit the first 2.5 s of data 
in Fig. [3] after fixing L'2 to its value scaled from 709 G 
and find if = (9.5 ± 0.5) x lO-^^cm^/s at 810 G. The 
slight disagreement in Li atom number at long times may 
be due to a small (< 10%) inequality in our spin mixture 
composition, which the model does not take into account. 

The qualitative features of both spectra in Fig.[T] can 
be theoretically reproduced by using field-dependent re- 
action coefficients scaled from our measured values at 709 
and 810 G. However, a full quantitative comparison will 
need to take into account the theoretical deviations from 
scaling behavior in the unitary regime as well as experi- 
mental variations in the initial temperature, and is open 
to future investigation. 

By extending the forced evaporative cooling step, lower 
temperature mixtures can be produced where bosonic 
^''^Yb shrinks to a size smaller than the Fermi diame- 
ter of the ^Li cloud. Such experiments at 834 G yield 
Tu/Tf ~ 0.25 with NYh = Nu = 2.5 x 10"^. Here, the 
estimated volume of the Yb sample is ~ 0.3 of the Li 
sample volume, compared to 3.3 in the classical regime. 
The mixture is thus also capable of achieving the oppo- 
site regime of a second species being immersed inside a 
strongly interacting quantum degenerate Fermi gas, sim- 



ilar to earlier studies in the K-Li mixture [7|. 

Our experiments with the Yb-Li mixture near a Fes- 
hbach resonance demonstrate effects of an additional 
species on chemical reaction rates in the microKelvin 
regime. We observe a long lifetime for Feshbach 
molecules, even in the absence of Pauli blocking. Our 
demonstrated stability of the mixture near the unitary 
regime of the resonance opens various possibilities of 
studying strongly interacting fermions immersed in a 
bath species or being interrogated by a small probe 
species. Future experimental opportunities include real- 
izations of non-equilibrium states, and studies of super- 
fluid properties, for instance by controlled relative motion 
between the two species. Finally, our results constitute 
an advance in the manipulation of ultracold mixtures of 
alkali and alkaline-earth-like atoms [igI . Isil . l32j . These 
mixtures may be used to produce quantum gases of para- 
magnetic polar molecules which are promising for diverse 
applications such as quantum simulation of lattice spin 
models [11], tests of fundamental symmetries [13], and 
probes of time variations in fundamental constants (35| . 
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